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Abstract 

An inverse design procedure was developed for the aesign 
of a mistuned rotor . The design requirements are that the 
stability margin of the eigenvalues of the aeroelastic system 
be greater than or equal to some minimum stability margin, and 
that the mass added to each blade be positive. The objective 
was to achieve these requirements with a minimal amount of 
mistuning. Hence, the problem was posed as a constrained 
optimization problem. The constrained minimization problem 
was solved by the technique of mathematical programming via 
augmented Gagrangians . The unconstrained minimization phase 
of this technique was solved by the variable metric method of 
Broyden, Fletcher, and Shanno. 

The bladed disk was modelled as being composed of a rigid 
disk mounted on a rigid shaft. Each of the blades were 
modelled with a single tosional degree of freedom. Adamcyzk 
and Goldstein's linearized aerodynamic model for the unsteady 
moment coefficients in a supersonic cascade was applied at the 
typical section. The resulting non-self-adjoint eigenvalue 
problem is of the form Aq = XBq. The eigenvalues and 
eigenvectors of this eigenvalue problem were found by a 
fourth-order Runge-Kutta line integration of the derivatives 
of the eigenvalues and eigenvectors. 

It was shown that mass mistuning does not introduce 
damping into the system, and that a necessary but not 
sufficient condition for stability is teat tne blade be self 
damped. The results of the optimization showed that an 
optimally mistuned rotor can achieve a given stability margin 


for a much lower level of mistuning than alternate mistuning. 
However/ it was shown that optimal mistuning is sensitive to 
errors in mistuning. Small errors in the implementation of 
optimal mistuning can severely reduce the gains in stability 
achieved by optimal mistuning. Alternate mistuning, on the 
other hand, is relatively insensitive to errors in mistune. 
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1 . INTRODUCTION 

The aaroelastic phenomenon of flutter in a transonic 
rotor is unusual in the field of airfoil aeroelasticity in 
chat the instability is not of the frequency coalescence type 
which is common in isolated airfoils. The instability in a 
rotor is due to the destabilizing effect of the cascade 
aerodynamics. The motion of a single blade will cause 
unsteady aerodynamic forces on all the other blades. These 
forces can cause the rotor to be unstable. It has been shown 
by several researchers that the effect of mis tuning, that is 
the altering of the natural frequencies of the blades of a 
rotor, generally has a benificial effect on the stability of 
the rotor [1-3]. In the past, however, the analysis of 
mistuned rotors has been limited to the determination of the 
aeroelastic behavior of a rotor whose mistuning has been 
specified. In this investigation, the inverse problem is 
solved. A method is presented for the determination of the 
mistuning arrangement which provides the greatest stability 
for the least amount of mistuning. 

Recently, Kaza and Kielb [1,2] have used a 
bending-torsion coupled model of the blades to examine the 
effects of raistuae on subsonic and supe 'sonic rotors . They 
showed that if the elastic axis is at the midchord of the 
blade, and the first bending frequency and the first torsional 
frequency of the blades are not nearly equal, then the flutter 
speed predicted by a single torsional mode agrees well with 



the results of the two degree of freedom bending-torsion 
model* This indicates that the flutter of rotors is not of 
the clasical bending-torsion frequency coalescence type. 
Furthermore, their work has demonstrated that the effect of 
raistuning is generally stabilizing. 

The phenomenon of flutter in transonic rotors is 
dependent on the unsteady aerodynamic forces acting on the 
blades of the rotor. Several authors have studied the 
unsteady forces acting on blades in supersonic flows [4-6]. 
The aerodynamic conditions modelled were steady uniform flows 
over flat, nonturning airfoils with unsteady perturbations of 
the flow due to the motion of the airfoils. The flow was 
modelled as irrotational , isentropic, and two-dimensional. 
The blades undergo simple harmonic motion but are phased such 
that the motion of the blades can be described in terras of 
travelling waves. This representation of the blade motion is 
useful for deriving these unsteady forces. However, it will 
be shown that this representation is not as useful for 
studying the mis tuned aeroelastic behavior of the rotor. 

The unsteady aerodynamic forces can, through a Fourier 
transformation, be converted into an influence coefficient 
form. The coefficients indicate the force felt on a reference 
blade for a unit amplitude displacement of any other blade on 
the rotor. In this form, one can clearly identify the origin 
of the forces acting on a blade. As early *3 1969, 
Samoylovich [7] used such transformations to determine the 


influence coefficients of an Infinite cascade. In 1950, 
Srinivap^.n [8] used a similar transformation to study the 
influence of mistuning on blade torsional flutter of a 
shrouded fan. In that same year, Hanamura, Tanaka, and 
Yamaguchi [9] used the inverse of this transform to convert 
experimental data from individual blade generalized 
coordinates to travelling wave coordinates. Recently, Kaza 
and Kielb [1,2] have used a similar transformation to 
transform the mass and stiffness matrices from the individual 
blade coordinates to the travelling wave coordinates to study 
the effects of raistune on the flutter and forced response of 
rotors. 

In Chapter 2 of this report, the basic theory of rotor 
aeroelasticity is reviewed. The blades are modelled with one 
degree of freedom per blade to study the effects of structural 
mistuning on the aeroelastic phenomenon of flutter . The 
mistuned rotor stability problem is cast as a matrix 
eigenvalue problem. By making use of the properties of the 
eigenvalue problem, it will be shown that structural mistuning 
does not introduce damping into the system, but rather makes 
use of the damping already present in the system to delay the 
onset of flutter. Hence, there are limits to the usefulness 
of mistuning as a mechanism for the preventition of flutter. 

In Chapter 3, a method of determining the eigenvalues and 
eigenvectors of the mistuned system is presented. The first 
step in the process is to determine the derivatives of the 
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eigenvalues with respect to the mistune parameters. Once 
these derivatives are known, one can then integrate these 
derivatives to another ra,%stuned state to determine the new 
eigenvalues and eigenvectors. This integration is carried out 
using a fourth-order Runge-Kutta algorithm. It was necessary 
to develop this method of evaluation of the eigenvalues for 
use with the inverse design procedure discussed in Chapter 4. 
Conventional methods of evaluating the eigenvalues of such a 
system (such as the inverse power method or the QR method) are 
shown to be unacceptable for use in the inverse design 
procedure . 

The inverse design procedure is presented in Chapter 4. 
When designing a rotor which will not flutter, one would like 
to minimize the amount of mistuning required to achieve a 
flutter free system. Por this reason, the inverse design 
procedure was posed as a constrained minimization problem. 
The minimized quantity is the amount of mistune in the rotor 
and the constraint insures an adequate stability margin. This 
inverse design problem can be divided into two parts; the 
definition of the problem statement; and the efficient 
solution for the constrained minimum. As part of the problem 
definition, a cost function must be defined which represents 
the level of mistune which is present in the rotor. The 
design specifications are then cast as constraints. For 
instance, one such design requirement is that the rotor be 
free of flutter at its aeroelastic design point. The second 


part of the problem ia to find the solution for this 
constrained minimization problem. The solution is defined as 
that mistune pattern which minimizes the cost function while 
satisfying ail the design requirements. The technique used to 
solve this problem is mathematical programming via augmented 
Lagrangians [10,11]. 

Finally, in Chapter 5, the results of the optimization 
are presented. It will be shown that optimal mistuning can 
achieve the design requirements at a relatively low level of 
mistuning compared to alternate mistuning. However, several 
other issues need to be addressed before the designer can 
implement these optimal mistune patterns- For example, the 
optimal mistune patterns found in this investigation are very 
sensitive to errors in implementation. If the rotor is not 
raistuned very precisely, the rotor will lose the stability 
margin gained in the optimal mistuning- This and facets of 
the problem of practical implementation will be discussed, 
leading to a realistic assessment of the practical value and 
realizable optimization of rotor mistuning. 
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2. AEROELASTIC EQUATIONS OP MOTION FOR A BLADED DISK 

In this chapter, the equations of motion which govern the 
aeroelastic behavior of a bladed disk are developed and the 
implications these equations have on stability are examined. 
The bladed disk is modelled as a rigid disk mounted on a rigid 
shaft. Each of the N flexible blades are aerodynamicaily 
identical but may have small differences in structural 
properties from one blade to the next. This structural 
nonuniformity is known as mistuning . Mistuning may be of a 
statistical nature due to manufacturing tolerances or material 
differences from blade to blade, or it may be designed into 
the rotor by introducing deliberate changes in blade materials 
or dimensions . 

2.1 STRUCTURAL DYNAMIC MODEL OF THE BLADED DISK 

The disk considered in this investigation is assumed to 
be rigid and mounted on a rigid shaft rotating at rotational 
speed Each blade is assumed to have a single torsional 
degree of freedom. In this investigation, this degree of 
freedom was taken to be the first torsional mode of the blade. 
The structural model does not include the effects of disk 
flexibility, rotation, or blade to blade coupling through 
shrouds. However, the principles of mistuning can be extended 
to include such effects as discussed by Kaza and Kielb £12]. 
In the absence of aerodynamic forces, each blade is uncoupled 
from every other blade and behaves like an uncoupled 
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oscillator. In this typical section analysis , the resulting 
equations of motion are 




o. rt-i 


( 2 . 1 ) 


where is the typical section moment of inertia, is 
the first natural torsional frequency of the ith blade, and 
q^ is the torsional degree of freedom of the ith blade about 
the elastic axis as shown in Figure 2.1. In the presence of 
aerodynamic forces, the blades are recoupled since the forces 
acting on one blade depend on the time. history of all the 
other blades. 


£ *j i + -i k , «7 i, 

i - 0, tjl, V-/ 


( 2 . 2 ) 


where M * the relative Mach number 

k = the reduced frequency of vibration 
a = the solidity of the rotor at the typical 
section 

£ = the stagger angle of the typical section 
a = the location of the pitch axis 


In matrix form this becomes 
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The right hand aide of Equation (2.3) includes the effect 
of both forces due to external forcing and forces arising from 
the motion of the blades. Or said another way, the blades are 
subjected to both inhomogeneous and homogeneous forces . The 
stability of this system is governed by the homogeneous terms. 
Therefore, when analyzing the aeroelastic stability of the 
system, only the motion dependent forces need to be included. 

To determine the solution to Equation (2.3), it is 
assumed that all the blades undergo simple harmonic motion. 
This assumption is made for two reasons. First, unsteady 
aerodynamics coefficients are derived assuming the airfoils of 
a cascade undergo uniform harmonic motion. Second, even 
though such an analysis does not strictly give the damping 
ratios of the different modes, such an analysis will correctly 
predict the neutral stability mode3 since the blades undergo 
simple harmonic motion at the stability boundary. With the 
assumption of simple harmonic motion, the displacement of the 
blades is expressed as 

Cj’ « t ) } l* °j *j i-j"- , N ~' (2.4) 

Substitution of Equation (2.4) into the equations of motion 
results in the eigenvalue problem 



(2.5) 



The aerodynamic influence coefficients, CL], are as yet 
undetermined. However, by the symmetry of the cascade, it 
must be that CL] is a circulant matrix C3]. That is, all the 
columns of [L] are identical except that each succesive column 
is shifted down by one row since the effect of blade 2 on 
blade 1 must be the same as the effect of blade 3 on blade 2, 
etc. Of course, [L] is still a function of w since the 
unsteady aerodynamics are a function of the reduced frequency 
of vibration, k = u)b/U. Dugundji and Bundas C13] have shown 
how to approximate the unsteady aerodynamic forces over a wide 
range of frequency k using Pade approximates and have included 
this effect into a standing wave flutter analysis. However, 
for the present analysis it is assumed that [L] is independent 
of u since, for small amounts of mistune, the reduced 
frequency of the various modes varies by a small amount. 


2-2 NONDIMENS IONAL FORM OF THE EIGENVALUE PROBLEM 

Some insight into the problem of flutter can be gained by 
examining the eigenvalue equations in nondimensional form. 
Toward this end. Equation (2.5) is divided by I R ui R 2 where & R 
is the first torsional frequency of the nominal reference 
blade and i_ is the typical section moment of inertia, 
the nondimensional eigenvalue problem becomes 


Then 
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r\ 1 

f-1 


N 

- 


= xi 1 


r N 






( 2 . 6 ) 


where e^=* the fractional mass raistune of the ith blade 
6 ^" the fractional stiffness raistune of the ith 
blade 

=• (u/to R “ the nondiraensional eigenfrequency 

U => the nondimensional mass of the blade compared 

to the mass of the surrounding fluid 
2 

* m/ffpb 

r * the nondimensional radius of gyration 
a ^j~ I/mb^ 

m = the mass of the typical section of the 
nominal blade 

In Equation (2.6), the matrix [L] appears as an apparent 

inertial term. This is a result of the choice of notation 

since the aerodynamic forces were assumed to be proportional 
2 

to w CL] . If instead the aerodynamic forces were chosen to 

2 

be proportional to [l] , where ui R is the natural frequency 
of the nominal blade, then [L] vrauld appear as an apparent 
stiffness term. However, both are good approximations at the 
stability boundary for small levels of mistune. Dugundji and 
Bundas C13] have used Bade approximates to obtain better 
approximations to the aerodynamic forces both close to and 
away from the stability boundary and for a large range of 
reduced frequency. 

Note that this eigenvalue problem is of the form 
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(2 ‘” 

This form of the eigenvalue equation is slightly more 
difficult to work with than the standard eigenvalue .form, 
Aq » Xq. Some useful properties of Equation (2.7) are 
developed in Section 3.1. 

In general, the eigenvalues X will be complex and hence, 
ft will also be complex. The fact that 0 ia not purely real 
but will have a small imaginary part violates the assumption 
of simple harmonic motion. Strictly speaking, these 
eigenvalues do not give the correct damping ratio of the 
different eigenmodes. However, the assumption of simple 
harmonic motion is only violated slightly since the imaginary 
part of ft is usually very small compared to the real part. 
Bundas [14] has shown that for small damping ratios, the 
damping ratio predicted by assumming simple harmonic motion 
gives a good approximation to the actual damping ratio. 

This approach of assuming simple harmonic motion, only to 
find that the eigenvalues do not represent simple harmonic 
motion, is similar to the traditional V-g diagram analysis 
often used in analyzing the stability of isolated airfoils. 
The damping ratio of each eigenvalue found from Equation (2.6) 
is interpreted as the damping ratio which must be subtracted 
from that mode so that the mode will undergo simple harmonic 
motion . 

An alternative formulation to the eigenvalue equation 


given in Section 2.2 is to pose the eigenvalue problem in 


terms of travelling wave coordinates. Formulation of the 
nondimensional eigenvalue problem in this fashion gives 

[eflVvJMW = 

( 28 ) 

where [E] is a linear transformation matrix (see Section 2.3), 
and q e is a travelling wave with an interblade phase angle 
of 3 n = 27rn/N. This form has been used by other investigators 
to study the effects of mistune on rotor stability and forced 
response [1,2]. The principle advantage to this form is that 
the eigenvalues and eigenvectors of the tuned system are 
readily determined. For and 6^ = 0, the eigenvalue 

problem becomes 

[ r ]{v] = < 2 - 9 ) 

The characteristic equation is then 

det(ci] ‘ ° (2-io) 

Since all the matrices in Equation (2.9) are diagonal, the 
determinant is equal to zero if and only if one or more of the 
diagonal entries is equal to zero. This is equivalent to 


k - 


= V, 3 -; 


<[’ + = 1 


N - / 


( 2 . 11 ) 
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Hence, the eigenfrequenci.es are 

( 2 . 12 ) 

It is seen from Equation (2.12) that the tuned rotor will 
flutter if any one of the unsteady moment coefficients l e has 
a positive imaginary part. This formulation of the eigenvalue 
problem is not as useful, when and <5^ are not equal to 
zero. In the mistuned case, the matrices £E] ^ [1+ e^]£E] and 
[E]” 1 £1+ 6^ J£E J will, in general, be fully populated. 

2.3 DETERMINATION OF THE AERODYNAMIC INFLUENCE COEFFICIENT 

MATRIX i 

i 

When describing simple harmonic motion of the blades, t> 

there are many equally valid generalized coordinate systems 

t 

one can use to represent their motion. The three common 
choices are the travelling wave representation, the standing \ 

I; 

wave representation, and the individual blade coordinates 
where each blade is represented by its own degree of freedom. ; 

Of course, one may easily transform from one coordinate system 
to another by simple coordinate transformations. The choice 
of coordinate system will generally depend on the nature of 
the problem being solved. For instance. When analyzing a 
flexible disk, the natural choice is the standing wave 
coordinate system since for the tuned rotor, the eigenmodes 
are standing waves which can be classified by the number of 
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nodal diameters and nodal circumferences . 

‘ When deriving unsteady aerodynamic forces which act on a 
blade, it is easiest to work in travelling wave coordinates. 
However, describing the force acting on a reference blade as a 
function of interblade phase angle does not yield a physical 
understanding of the origin of these forces. One would like 
to look at the unsteady air loads as influence coefficients, 
i.e., forces which act on a given blade due to the motion of 
other individual blades. Looking at the forces in this manner 
gives a clearer picture of the mechanisms of flutter. 

Consider the cascade shown in Figure 2.2. Suppose there 
is a travelling wave moving through the cascade with an 
interblade phase angle Of 2"rrn/N. The deflection of the 
ith blade is then 


V = % 




(2.13) 


As has been defined, the travelling wave moves in the 
direction of rotor rotation with a wave speed of tu/n. Of 
course, it is equally valid to consider this motion as a 
backward travelling wave with an interblade phase angle of 
2tt - To completely describe all possible combinations of 

sinusoidal motion of the N individual blades requires N 
interblade phase angles, 


ft*o 


(2.14) 


bi 
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This is conveniently expressed in matrix form as 


M ■ l‘]W*'"‘ 


( 2 . 15 ) 


where the matrix [E] is 


u= 




9 * • 


^*l 4 » ^>,1 4 . . 

^'Si • • 


* * 


L 




• • 


6 ***V»I 


( 2 . 16 ) 


where 


u ■ 


The unsteady aerodynamic moment acting on the zeroth 
blade for a given travelling wave with amplitude and 
interblade phase angle 8 n is 


n 




jal't 


( 2 . 17 ) 


But the force acting on the ith blade due to travelling wave 
mode n is just the force acting on the zeroth blade but 
shifted in phase by 8 n i . Summing the contributions made from 
all the interblade phase angles gives 

( 2 . 18 ) 

n*o 



t 

r 
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This is expressed more succinctly in matrix form. 




(2.19) 


Next, Equation (2.15) is used to eliminate qg from Equation 


n 


(2.19). This gives the desired result of describing the 
unsteady aerodynamics in terms of the motion of the individual 
blades 


M - ^e^][X][<f (<?<} 


( 2 . 20 ) 


where the entries of the matrix [E] ^ are given by 


, - }*b£1 

f~ • J- e * 

N 


( 2 . 21 ) 


So that finally the influence coefficient matrix CL] is given 
by 


H-MvW 


( 2 . 22 ) 


Multiplication of the right hand side yields that CL] is 
in fact a circulant matrix. 




L# i-v-i L »/- «■ * * • L , 

^ 1 • L// -i .*• L t 

L < L .... 


• • 


,L-a/-i Ln-x « • • L * 


(2.23) 
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where 


v - ( 


^ = ~kl2 h* £ 


rt»o 


j jJT-Xfl. 

N 


(2.24) 


Equation (2.24) shows that L K is just the Kth 
coefficient in a discrete Fourier series representation of 

x e * 

n 

a -jULKn. . 

'1Z L « e - ~ (2 - 2s > 

K»o 


So for example, if [L] is a tridiagonal matrix, then the 

physical interpretation is that only the two blades adjacent 

to a given blade and the blade itself have any direct effect 

on the blade. Then from Equation (2.25), l fi will contain 

only cos(0), cos(0), and sin(B) components. In other words, 

the plot of the real and imaginary parts of lg when plotted 

as a function of 0 will have a D.C. offset and will have 

components of cos(0) and sin(0). If the influence coefficient 

matrix has a larger bandwidth, then l e will contain higher 

p n 

harmonics of cos(0) and sin{6). 

As early as 1969, Saraoylovich C7] used a similar Fourier 
transform to convert from travelling wave coordinates to 
individual blade coordinates. In 1980, Hanamura, Tananka, and 
Yamaguchi C9] reversed the process. They experimentally 
measured the influence coefficients of a cascade of blades in 
incompressible flow by vibrating a single blade and measuring 
the resulting forces on all the other blades. They then 
converted these forces to interblade phase angle coordinates 


by the use of the Fourier series representation. Equation 
(2.25) 


2.4 IMPORTANCE OF BLADE SELF DAMPING 

As stated earlier, the matrix [L] is circulant. Hence, 
every term on the diagonal is equal. This terra, denoted by 
Lq , reflects the aerodynamic effect the motion of the blade 
has on itself* To show the importance of this term, consider 
the case of mass mistune only (i.e., no stiffness mistune, 
<5^ * 0) . Then the eigenvalue problem can be written as 
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(2.26) 


Here, 1 /ll represents the eigenvalues of the matrix on the 
right hand side of Equation (2.26). Next, making use of the 
familiar matrix property that the sum of the eigenvalues of a 
matrix is equal to the trace of the matrix, it must be that 
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(2.27) 


Recall that the nominal reference blade in the absence of 
aerodynamic forces vibrates at the nondimensional frequency 
= l. Since the unsteady aerodynamic forces are small 
compared to the elastic and inertial forces, we expect that 
the complex nondimensional eigenfrequencies fi will be very 
close to unity. For convenience, let = s. This gives the 
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familar s-plane interpretation of the poles or 
e igenfrequencies . If any pole lies in the right half of the 
complex plane, then the system is unstable. 

In the remainder of this section, necessary conditions 
for stability will be derived. For stability, all poles must 
lie in the left half of the complex s-plane. Since ft is 
approximately equal to unity, s will be nearly equal to j. 
Let s be represented by s * j + n where s is a complex number 
much less than unity. Then 


n 


= 'J* 


(2.28) 


so that 


i _ - 1 

V 

= / ' Ajs * 0CS‘) 


(2.29} 


Next the centroid of the poles in the complex plane, denoted 
by <s>, is found. Substitution of Equation (2.29) into (2.27) 
gives 


£2>**iV 

k 




< 3 > 25 



(2.30) 


The centroid of the poles is <s> so that from Equation (2.30) 


R*<s> “ir Em (^v) 

r "< J > » t-i 

Equation (2.31) shows, at least for small amounts of 
mistune and large ur , that the centroid of the poles lie in 
the left half of the complex plane if and only if Im(Lg) is 
less than zero. This is equivalent to saying that a necessary 
but not sufficient condition for stability is that the blades 
must be self damped. Or said another way, if all the blades 
were perfectly rigid except for a single flexible reference 
blade, that blade must not flutter. This condition is the 
fundamental limitation to the usefulness of mistuning as a 
mechanism for stabilizing a bladed disk. Note that Re<s> is 
independent of e^. Physically, this is because mistuning 
does not introduce any damping into the system. 

On the other hand, Im<s> does depend on The 
lowering of the centroid by < e >/2 in the complex plane just 
reflects the fact that the natural frequency of each blade is 
inversely proportional to the square root of 1 + e i . 
Similarly, the term Re(I^/ur^) can be thought of as an 
effective mass added to the blade due to the unsteady 
aerodynamics, or alternatively, as an aerodynamic destiff ening 
of the blade. 


2.5 EXAMPLES OF TUNED AND MISTUNED ROTORS 

As an example, consider the unsteady loads acting on a 

typical section of a blade. The aerodynamic model usee in 

this example and throughout this investigation is that of 

Adamczyk and Golds tain's [4]. In this model, linearized 

theory is used to obtain the unsteady loads on flat plate 

airfoils undergoing small amplitude harmonic oscillations. 

Shocks are modelled as isentropic Mach waves and there is no 

steady pressure rise across the cascade. The typical section 

has a single torsional degree of freedom and pitches about its 

midchord. In this example, the reduced frequency, k, is equal 

to 0.495, the Mach number, M, is equal to 1.317, the solidity, 

a, is equal to 1.409, and the number of blades, N, is equal to 

14. These rotor parameters are tabulated in Table 2.1. The 

unsteady moments for this model are plotted in Figure 2.3 as a 

function of interblade phase angle. Note that for interblade 

phase angles 25.71°, 51.43°, 77.14°, and 102.86°, the 

imaginary part of 1 0 is positive. Therefore, the rotor will 

flutter at this reduced velocity and Mach number in its tuned 

position. However, the average value of 1 g is less than 

p n 

zero. This indicates that the centroid of eigenvalues lies in 
the left half of the s-plane and hence, it may be possible to 
achieve aeroelastic stability through the use of mistuning. 

Next, Equation (2.24) is used to transform the unsteady 
moments from interblade phase coordinates to influence 

coefficient form. These influence coefficients are plotted in 



Figure 2.4. Several interesting features of the influence 

coefficients deserve mention. First of all, the term has 

a negative imaginary part. This again reflects the fact that 

the average value of the imaginary part of 1 0 is negative, or 

n 

that the blade self damping is stabilizing, and hence, the 

potential for stability exists. Secondly, the largest 

coefficients are seen to be Lg, L^, and 1 ^ 3 . It is clear that 

only the near field neighbors of a given blade exert a 

significant influence on the blade. The and L^terms are 

the first off diagonal terms of [Lj. Examination of Equation 

(2.25) reveals that a dominantly tridiagonal [L] results in an 

lo which has a strong first harmonic dependence in S • 
p n n 

Figure (2.3) shows that this is in fact the case. 

The previous discussion suggests that alternate mistiming 
may be an effective mistuning arrangement since one would 
expect such a pattern to reduce the influence adjacent blades 
have on each other. Figure 2.5 shows the eigenvalues of the 
tuned system plotted in the s-plane. Figure 2.6 shows the 
eigenvalues of the alternately mass mistuned rotor. The even 
numbered blades have e^= 0 while the odd numbered blades have 
e i =» 0.1 . This mistune pattern does in fact stabilize this 

rotor. However, it remains to be seen whether or not 
alternate mistuning is a near optimal mistuning arrangement. 
This discussion is deferred to Chapter 5. 
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3. EVALUATION OP EIGENVALUES AND THEIR DERIVATIVES 

In the previous chapter, the stability of a rotor was 
found to be dependent on the eigenvalues of the aeroelastic 
eigenvalue problem. In the next chapter, an optimization 
procedure will be described. The goal of this optimization 
will be to find the mistune pattern which provides the most 
stability for the least amount of mistuning. In this 
optimizaton, it will be necessary to identify each of the 
eigenvalues uniquely. The technique used is this 
investigation to evaluate the eigenvalues while insuring 
unique identification will now be described. 

Consider for the moment the tuned rotor under the 
influence of unsteady aerodynamic forces. For this system, 
there are N eigenvalues. Corresponding to each eigenvalue is 
an eigenmode. The nth eigenmode is a travelling wave with an 
interblade phase angle of 2im/N. However, as soon as the 
system is mistuned, the eigenmodes are no longer pure 
travelling waves of a single interblade phase angle. The 
problem is to in some sense identify each eigenvalue of the 
subsequent mistuned system with the eigenvalues of the 
original tuned system. 

As an analogy, consider that the mistune introduced into 
the system is a gain and that the eigenvalues in the complex 
s-plane are the poles of a control system. As the gain is 
increased, i.e., as the level of mistune is increased in a 
continuous way, the poles, or eigenvalues, trace out root loci 




in the complex plane. Hence, in some sense, one of the 
eigenvalues of the mistuned system will be identified with the 
nth eigenvalue of the tuned system since they both are on the 
same root locus. 

It is not, however, an easy task to identify in this 
sense the nth mistuned eigenvalue. In general, the eigenmode 
of the nth mistuned eigenvalue will be composed of travelling 
waves of all N interblade phase angles. It is not generally 
possible to identify the the loci of roots that a particular 
root belongs by inspection of the eigenmodes. For these 
reasons, routines such as EISPACK [15] are not acceptable for 
the purposes of this investigation. 

Therefore, before the optimization of Chapter 4 could be 
performed, a scheme had to be developed which could evaluate 
the eigenvalues of the aeroelastic eigenvalue problem while 
uniquely identifying the root locus to which each eigenvalue 
belongs. That method of determining the eigenvalues and 
eigenvectors of the aeroelastic equations of motion is the 
topic of this chapter. In Section 3.1, it is shown how the 
derivatives of the eigenvalues and eigenvectors with respect 
to the mistune parameters are determined. Once the 
derivatives of the eigenvalues and eigenvectors are known, one 
can integrate these quantities from one mistune level to 
another through the use of an explicit fourth-order 
Runge-Kutta operator. The details of such an integration are 
given in Section 3.2. Finally, an illustrative example is 
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presented in Section 3 . 3 which demonstrates the use of this 
technique . 


3.1 NONDEGENERATE PERTURBATION THEORY APPLIED TO THE 
NON-SELF-ADJOINT EIGENVALUE PROBLEM 

In this section, the derivatives of the N eigenvalues of 
the aeroelastic eigenvalue will be found with respect to the 
mis tune parameters . First, some properties of the 

eigenvalue problem will be reviewed. These properties will 
then be used in a perturbation analysis of the eigenvalue 
problem. At a given level of mistune, a pertubation parameter 
will be introduced into the problem. The end result of the 
perturbation analysis will be that the derivatives of the 
eigenvalues and eigenvectors of the mistuned system will be 
“known with respect to the single perturbation parameter . 
Finally, the derivatives will be generalized to produce the 
derivatives of the eigenvalues and eigenvectors with respect 
to each and every mistune parameter. 

The eigenvalue problem 

MM " (3.D 

is known as a non-self-adjoint eigenvalue problem whenever one 
or both of the matrices A and B are not symmetric. Under 
these circumstances, the eigenvectors of Equation (3.1) are 
not in general the same as the eigenvectors of the adjoint 
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eigenvalue problem [16] i 

L ~ °lim 1 [b] (3.2) 

However, both eigenvalue problems have the. same set of 
eigenvalues since they both have the same characteristic 
equation 

(LA] - ^ £ aj) = o (3.3) 

The eigenvectors of Equations (3.1) and (3.2) are known as the 

right and left eigenvectors respectively. To obtain the 

relationship between the right and left eigenvectors. Equation 

T 

(3.1) is premultiplied by the mth left eigenvector q^ , 

'Vf 

Equation (3.2) is postmultiplied by the nth right eigenvector 
T 

q Rn • and the resulting equations are then subtracted one from 

<v * 

another to obtain 

(*hn\ ~ ® (3.4) 

It is assumed at this time that there are no repeated 
eigenvalues. Therefore, the only way for Equation (3.4) to be 
nontrivial ly satisfied is if 


l cj tm j [a]|c[^j - 


o 


m ¥■ n 


(3.5) 


This is known as the biorthogonality condition I*. 16] and will 
be useful in the perturbation analysis which follows. 

In Chapter 4, a procedure for optimally mistuning a rotor 
is described. As will be shown/ it will be necessary to 
evaluate the derivatives of the eigenvalues and eigenvectors 
to perform this optimization. The remainder of this section 
deals with the determination of these derivatives through a 
perturbation analysis. 

Suppose that the eigenvalues and eigenvectors of the 
general eigenvalue problem (3.1) are known. Next the matrices 

A and B are changed slightly. These changes are due to the 

'V 'b 

introduction of- a perturbation parameter y. It is assumed 
that the matrices A and B can be expressed in terms of a 

^ 'b 

Taylor series in y . Then 

4 - A m + TA*' + rW"> ••• 

a - b ” «. fB 0> +■ r’s w f 

where A^ , B^, A^ , B^, and so on are known quantities. 

'b 'b 'b 'b 

Next it is assumed that the eigenvalues and eigenvectors of 
the problem can also be represented in terras of a Taylor 
series in y ; 


(3.6) 

(3.7) 
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V-/ 






(3.10) 


where all the terms appearing in Equations (3.6) to (3.10) are 
of order unity except for, Y. The terms like a i n c TRi^ in 
Equations (3.9) and (3.10) represent a convenient change of 
basis for representing the changes in the eigenvectors. 
Substitution of Equations (3.6) to (3.9) into Equation (3.1) 
gives the following asymptotic relationship: 


(A m \ f 6"* ••-)(& + ~ 


Collecting terms in Equation (3.11) of equal order gives 


*<»>*> 0) _**) m) J/ .»> <•> 

5 p* ”* *<• s '*'T(~a 

+ 5 ^ S p. * V § T,°-J*pt) 

+ o(r x ) 


(3.12) 


For this asymptotic relationship to hold true, each 
coefficient of every order of y must vanish. This results in 
the equalities 
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Equations (3.13) and (3.14) give the zeroth and first 
order perturb - ' ' ..on relationships between the known quantities 
A (0)^ , B^ and the unknown quantities , 

r\j >\j r\j i\j n 

X n^ ' q Rn*' a i n * Notice that the zeroth order equation is 
just the unperturbed eigenvalue problem, and A^®^ and q Rj ^ 

'V 

can be found from its solution. This indicates that this is a 
regular perturbation problem. 

To determine the first order changes in the eigenvalues. 
Equation (3.14) is premultiplied by and the 

biorthogonality condition is imposed to eliminate terms equal 
to zero with the result that 
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(3.15) 


Solving for the unknown quantity A^ 1 ^ gives 
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(3-16) 


Hence, if A and B are perturbed about some point A^ and 

r\, r\j *V> 

(0) 

B , the value of the nth eigenvalue can now be estimated 
^ 2 

with errors of order y by substitution of equation (3.16) 
into Equation (3.8). 

The perturbed eigenvectors are found in a similar fashion. 

/ a ) T 

Equation (3-13) is premultiplied by q' ' • Again, the 

v jp 

biorthogonality relations are used with the result that 
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<* JT „ . A <•) V 0) miT <•* . V*' w r «•»<*> 


/• l ( , ) r ,.»<.) (3.17) 



(3.18) 


The left eigenvalue perturbation terms are found in a 
completely analogous fashion. 



l“ V’ - l- ?’)£, 

oT-ipfjrG. 


(3.19) 


Note that Equation (3-18) is valid only if p ^ n and 

yi X^°l analysis is therefore known as nondegenerate 

perturbation theory [17] since the eigenvalues must be 

distinct. 

The scalar Sp n indicates the degree to which the pth 

eigenvector is coupled to the nth eigenvector through the 

introduction of a small perturbation. Note that eigenvectors 

with closely spaced eigenvalues are more easily coupled than 

( 0 ) ( 0 ) 

those whose eigenvalues are far apart due to the (X n -A ) 
term appearing in the denominator of Equation (3.18). From 
Equation (3.16) it is observed that the eigenvalue 
perturbations do not exibit this amplification due to the 
proximity of other eigenvalues. 

The range of validity of this perturbation analysis is 
that range for which yA^ 1 ^ * yq^ ' and Yq Ln 


are less than 
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order unity. Hence, from Equations (3.18) and (3.19), it is 

clear that for this to hold, it must be that Y is small 

compared to min(X^ 0 ^- X^ 0 ^), the spacing of the closest 

eigenvalues. It is this fact which ultimately limits the 

usefulness of such an approximation. 

An interesting feature of this first order perturbation 

analysis is that and b„„ are not determined. This is 

* nn nn 

due to the fact that there is some degree of freedom in 
choosing an eigenvector. Recall an eigenvector specifies a 
direction in N-space but not a length. Furthermore, to order 
y , a small perturbation of the eigenvector in the direction 
of the eigenvector only produces a change in direction of 


order y 4 


Hence, one is free to choose any values for a 


nn 


and b „ so long as they are no greater than 0(1). For 
nn - 

convenience, the values of a„„ and b„„ are taken to be zero. 

nn nn 


( 1 ) 


( 1 ) 


Although the first order coefficients X^ , q Rn ' and 

(1) ^ 
q Ln were derived from a perturbation analysis, Plaut and 

r \i 

Huseyin [18] have shown that these are in fact the first 
derivatives of the eigenvalues and the right and left 
eigenvectors with respect to the variable Y . In the mass 
mi3tuning problem, there are N independent variables denoted 
by the mass mistune vector e . The results of the previous 
section can be generalized to give the derivatives with 
respect to every mistune variable; 


* L s 


(3.20) 
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(3.21) 

%** 

(3.22) 


3.2 INTEGRATION OF EIGENVALUES AND EIGENVECTORS FROM THEIR 
DERIVATIVES 

In the previous section, the derivatives of the 
eigenvalues and eigenvectors were found with respect to the 
mistune parameters. Equations (3.20) through (3.22) make up a 
system of N + 2N first order, coupled, partial differential 
equations for the N unknown eigenvalues and 2N unknown 
eigenvectors. There are N independent variables e^. If the 
derivatives of the eigenvalues and eigenvectors are known 
everywhere along a line, then a line integration can be 
performed to evaluate the eigenvalues and eigenvectors at a 
point along this line. This method of solving for the 
eigenvalues and eigenvectors has two advantages. First, if 
the matrices A and B are changed slightly, then the new 

'V 'V- 

eigenvalues and eigenvectors can be evaluated by integrating 
over a short distance with less computational effort than the 
effort required to completely re-solve the eigenvalue problem. 
Secondly, the eigenvalues are automatically kept track of in 
the root locus sense. The problem of eigenvalue identity 
discussed in the introduction to this chapter is completely 
eliminated. This was the main reason for choosing this method 
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of eigenvalue evaluation. 

Now that the derivatives are known# they may be used in a 
line integration of the eigenvalues and eigenvectors. Assume 
the line of integration to be the line in parameter space that 
connects some old value of the independent variable to 

'V 

the new value e K+ ^ • The path of integration can then 
described in terms of a single scalar variable hi 



"(H.: W3 


(3.23) 


= H * k H K 

Where h is a scalar value which varies from zero to one along 

the path of integration. Equations (3.20) through (3.22) can 

2 

be reduced to a system of N + 2N coupled# first order# 
ordinary differential equations by the chain rule. 


d _ ol 6 , + sL^- +• 
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(3.24) 


The derivatives of the dependent variables along the line of 
integration and the distance along the path of integration can 
now be expressed in terms of a single scalar variable h. For 
example, the derivative of with respect to h is 
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where the terras 9X n /3e^ are given by Equation (3.20). 

The numerical scheme used to compute the line integrals 
is a fourth order Runge-Kutta scheme [19]. Truncation errors 
are then of the order of step size to the fifth power. This 
type of scheme is one of a class of explicit 
predictor-corrector methods- The method is carried out in 
four steps: 
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To integrate over larger distances, say from e. 


to 
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e II ' the P ath of integration is divided into some number of 
intervals, I, as shown in Figure 3*1. These intervals are 
picked sufficiently small to avoid significant truncation 
error but large enough to keep the amount of computation 
required at a reasonable level. The line integration is then 
carried out using I fourth order Runge-Kutta integration 
steps . 

This method of evaluating the eigenvalues and 
eigenvectors of the equations of motion is not efficient when 
the number of integration intervals is large. However, when I 
is small, this method requires computation time on the same 
order as EISPACK. Generally, when evaluating the eigenvalues 
and eigenvectors of the mistuned system in the constrained 
optimization procedure discussed in Chapter 4, 1 is less than 
or equal to two. Under these circumstances, this procedure is 
reasonably efficient. But more importantly, the identities of 
the eigenvalues are determined in the process of the 
integration. This is a requirement for the constrained 
optimization of the flutter problem- 

3 • 3 EVALUATION OF THE ALTERNATELY MISTUNED ROTOR BY 
RUNGE-KUTTA INTEGRATION 

In the previous section, it was shown how to evaluate the 
eigenvalues and eigenvectors of the aeroelastic eigenvalue 
problem by a Runge-Kutta integration. In this action, a 
numerical example of that technique is presented which 
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demonstrates the accuracy of this technique. 

The example will proceed as follows: first the exact 
eigenvalues of an alternately mistuned rotor will be found. 
Then the Runge-Kutta integration technique will be used with 
the path of integration divided into 1, 2, 3, and 10 intervals 
as described in the previous section. 

The aeroelastic properties of this rotor are listed in 
Table 2.1. This is the same 14-bladed rotor first introduced 
in Chapter 2. The rotor was alternately mistuned in mass. 
The masses of the even numbered blades was increased by 10 
percent while the roast: ?s of the odd numbered blades were 
unchanged. EISPACK was used to determine the exact values of 
the 14 eigenvalues. 

Next, the Runge-Kutta scheme was used to evaluate the 
eigenvalues of the system described above. The tuned 
eigenvalues and eigenvectors were use as the starting point of 
the integration since the eigenvalues and eigenvectors are 
easily determined at this point. Of coui 'e the accuracy of 
the integration will be determined by the step size of the 
intervals (inversely proportional to the number of intervals) 
used in the integration- The step size controls the level of 
truncation error. Hence, a smaller step size produces smaller 
total integrated errors. However# one must balance accuracy 
with computational efficiency. Integrations with small step 
sizes generally imply many steps and# hence, a large amount of 
computation . 
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Figures 3.2 through 3.5 show graphically the effect of 
step size on integration accuracy- Figures 3-2, 3.3, and 3.4 
show the results of the integration where the path of 
integration has been divided into one, two, and three 
intervals respectively. Notice that some of the eigenvalues 
evaluated by this approximate integration are very different 
than the exact eigenvalues as determined by EISPACK. This is 
due to the large step size used in the numerical integration. 
If the number of intervals is increased to 10, the approximate 
eigenvalues found from integration are virtually 
indistinguishable from the exact eigenvalues as indicated in 
Figure 3.5. 

The most useful feature of this method of evaluating the 
eigenvalues is not its efficiency since in fact for an 
integration which requires many intervals, this method is 
computationally very expensive. Its utility is in its ability 
to to retain the identities of the eigenvalues no matter how 
severe the mistuning. For this reason, this is the method 
used for evaluating the eigenvalues in the optimization 
procedure described in the next chapter . 


4. OPTIMAL MISTUNING OF A BLADED DISK TO PREVENT FLUTTER 

When miatuning a rotor to prevent flutter, one would 
like, in some sense/ to minimize the amount of mistuning while 
maximizing the stability of the rotor. Not only does such an 
optimization make the production of a mic tuned .totor more 
attractive but the results of this minimization may provide 
some insight into the nature of mistuning. 

In this chapter, such an optimization procedure is 
outlined. The objective of the optimization is to determine 
that mistune pattern which delivers the most stability for the 
least amount of mistuning. A cost function is devised which 
is a measure of the level of mistune in the rotor. This cost 
function is to be minimized subject to certain physical 
constraints on the problem. Although the desired effect is to 
maximize the stability of the rotor, it is more appropriate to 
modify the above optimization statement to make stability a 
constraint rather than a maximized quantity. The problem is 
then to minimize the level of mistuning in the rotor while 
meeting some minimum requirement of stability. Although these 
conditions may be stated very easily in mathematical terms, 
the solution to the problem is a computationally difficult 
task . 

Other authors have shown the useful benefits of alternate 
mistuning [1-3]. It has been suggested that the alternate 
mi3tune pattern may be nearly optimal C3]. The results of the 
optimization procedure presented in this chapter will be 
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presented in Chapter 5. It will be shown that in fact, 
alternate mistune is not nearly optimal, but does have other 
favorable properties that optimal mistuning does not always 
have . 

4.1 FORMULATION OF THE MISTUNING PROBLEM AS A NONLINEAR 
CONSTRAINED OPTIMIZATION 

In this section, the mistuning problem is formulated as a 
nonlinear constrained minimization problem. The steps to be 
taken are : first, chose an appropriate cost function which 
represents the level of mistuning in the rotor. Second, 
define the set of minimum requirements which must be met. 
These are known as constraints. Together, the objective cost 
function and the constraints form the constrained optimization 
problem. 

The first task in formulating the formal constrained 
optimization problem is to define the objective cost function. 
The objective cost function will be chosen to be a measure of 
the amount of mistuning in the rotor- As a simple case, 
suppose the cost of mistuning, <p , is chosen to be the 
absolute value of the mass fraction added to the blade with 
the greatest amount of mistuning. Then 

<pL I , J ••• i I (4.1) 


This cost function, although conceptually meaningful, is not 
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amenable to many nonlinear programming methods due to the 
discontinuities in the gradient of the cost at various 
locations in the domain* Note, however, that this cost 
function is equivalent to the cost function 



(4.2) 


where I is a positive even integer. For large but finite I 
this function closely approximates the cost function in 
Equation (4.1) but has no discontinuities anywhere. For this 
investigation, I was taken to be 4. 



(4.3) 


Originally, I was chosen to be 2, making the cost the 
root mean square value of the mass mistuning. In this case, 
the cost function is then just proportional to the length of 
the mass mistune vector e . However, it was found that this 
cost function, when optimized, does not penalize strongly 
enough large amounts of mistune in a single blade. The result 
is that the mass of one of the blades may become much too 
large to be practical. For this reason, I was increased to 4. 
This gave a more physically realizable distribution of 
mistuning. 

As will be shown, the optimization procedure outlined in 
this chapter requires that the derivatives of the cost 
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function and the constraints be evaluated. The gradient of 
this cost function is 



Having defined the cost function, the next step in 
formulating the constrained optimization problem is to define 
the appropriate set of constraints. As previously mentioned, 
the cost function will be minimized subject to minimum 
stability requirements. These stability requirements are 
interpreted as constraints. The measure of stability will be 
the damping ratio of the eigenvalues of the eigenvalue problem 
given by Equation (2.6). Recall that these damping ratios are 
not the physical damping ratios of the eigenmodes but rather 
are damping ratios in the sense of a V-g analysis. If the 
eigenvalue of the ith mode of Equation (2.6) is 
j^i 33 a i = u i + ^ v i' t * 1 ® 11 tb® damping ratio of the ith mode is 
given by 


1i = 


-14 




= Oj'jZj *-/ 


(4.5) 


One set of constraints is then that each and every mode of the 
rotor must have a damping ratio greater than or equal to some 
minimum damping ratio, 5 , as shown in Figure 4.1. This is 
expressed in the standard inequality form 
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9 1 L 6) - Jr t - - j ’ > O J l « Oj t , a, .. v V-/ 1 * 1 

which gives N constraints to be satisfied, represented by the 
functions 0^. 

Another set of constraints used in this investigation is 
that mass can be added to a blade to mistune the blade but 
mass cannot be removed.. This constraint was used because it 
was felt that a practical way to mistune blades would be to 
add masu to the tip of the blade or to reduce the stiffness of 
the root. Both of these methods tend to lower the natural 
frequency of the blade, in this problem, the stiffness of 
each blade was held fixed. It is believed that the natural 
frequencies of the individual blades dominates the mistuning 
effects. Hence, one may work with either mass or stiffness 
mistunings or a combination of the two with little difference 
in the resulting natural frequencies of the mistuned blades. 

This second set of constraints is expressed simply as 

= - 0 ) ‘ s Oj';!,.../-/ (4.7) 

which gives an additional N constraints to be satisfied. 

Together these two sets of constraints provide a total of 
2N inequality constraints for the minimization of a cost 
function of N variables. There were no equality constraints 
used in this investigation. 

The first derivative of the eigenvalue constraints 
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(Equation (4.6)) are 


Mi = f V,-iw - 

Uj (Ui^v i '? / * L >ij >tj J 
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(4.8) 


where ~ = - — '■ + Vl 3l2.*£ A 

where ^ z«/.'*rS)\ >4; » ) 

Wc - - / J u: l^tm: - 3X»« ; \ 

Z«*v;V > 4 * 6 i TiJ ' 

2 

where recall that X =* G and the derivatives of X are given 
by Equation (3. 20). The derivatives of the mass inequality 
constraints (Equation (4.7)) are 




(4.9) 


where is the Kronecker delta. 

To summarize, the purpose or the optimization procedure 
is to provide the most stability for a given rotor at its 
aeroelastic operating point for the lowest level of mistune. 
The formal mathematical statement of the problem is slightly 
different. The mistuned rotor must meet damping ratio 
requirements at a minimum cost which reflects the amount the 
rotor is mistuned. Furthermore, mass changes or mistunings in 
each blade must, in this investigation, be positive. 


4 2 FORMULATION OF THE CONSTRAINED OPTIMIZATION PROBLEM USING 

LAGRANGE MULTIPLIERS 


The general problem of nonlinear programming is to find 
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the minimum of some cost function subject to constraints. The 

cost function <t> is a function of the n variables x ; 

% 

$ * f (*) (4.10) 

In general, there may be two types of constraints which 
must be satisfied: equality and inequality constraints. 

Equality constraints are of the form 

?;(*) - ° j * tj ij ...j f> /> < h (4.11) 

While inequality constraints are of the form 

Oj(V^ 0 j j* (.V'O ? 7 < h (4.12) 

As indicated above, there may be any number of inequality 
constraints but there must :>e fewer equality constraints than 
the number of independent variables, n. The region in x-space 
where these constraints are satisfied is called the feasible 
region- The cost function <{> is said to have a constrained 
minimum [10] at x if there exists some positive a such that 

<j> ( l) £ $(%+*%) 

for all Ax in the set * 


(4.13) 
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Or said another way, a constrained minimum exists at a point x 
if x is in the feasible region and a small change in x in any 
direction in the feasible region causes an increase in If 

at such a minimum an inequality constraint function satisfies 
= 0, then that constraint is said to be active. If on the 
other hand, 9^ > 0, then the constraint is said to be 
inactive. 

Mathematically, it is convenient to add p+q additional 
variables to the problem by introducing Lagrange multipliers. 
Let the Lagrangian cost be 

^ -by r f j. yT q (4.14) 


Where 


v a the vector of equality Lagrange multipliers 

n the vector of inequality Lagrange multipliers 

It can be shown C10], except for rare circumstances known as 
abnormal cases, that the necessary conditions for a 
constrained minimum are given by 


VT = o 
Vi * ° > ‘ * 

Vj e r-° ] 
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j 


"j 


1 


(4.15) 

(4.16) 

(4.17) 

(4.18) 

(4.19) 


These conditions are known as the Kuhn-Tucker conditions. 

The addition of the product of the Lagrange multipliers 
and the constraint functions to the objective cost function <f> 
results in the Lagrangian cost J which has a slope of zero at 
the constrained minimum. In other words, the Lagrangian cost 
has a stationary point at the constrained minimum. Since the 
value of the constraint is zero at the constrained minimum, 
the Lagrangian cost J equals the cost <J> at that point. 
However, the stationary point which is a constrained minimum 
will not necessarily be a minimum of the Lagrangian cost. In 
fact, this point may be a saddle point or even a maximum of 
the Lagrangian cost. In other words, the Hessian of the 
Lagrangian cost {the Hessian is the matrix of second 
derivatives) will not necessarily be positive definite at the 
constrained minimum. This means that one cannot look for a 
stationary point by searching for a minimum of J. This 
problem is easily remedied by adding penalty functions to the 
Lagrangian cost, J, to form the augmented Lagrangian cost, 
J & . For the moment, consider only equality constraints. We 
wish the augmented Lagrangian cost to equal the Lagrangian 
cost everywhere along the constraints . But the Hessian of the 
new cost should be positive definite at the stationary point. 
This suggests adding a quadratic-like quantity to the 
Lagrangian cost which is zero everywhere along the constraint. 


Hence, we let 
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i * ft * 



(4.20) 


Where P must be greater than some minimum value to ensure 
that the Hessian of J a will be positive definite at the 
minimum. 

Figure 4.2 shows graphically the concept© discussed 
above. Consider the one-dimensional problem where 
<P (x) » 2 - x and ^(x) = x - 1. The cost, Lagrangian cost, 
and augmented Lagrangian cost functions are plotted versus x. 
By inspection the solution is x 1 and the Lagrange 
multiplier is easily found to be v =* 2. The Lagrangian cost 
is seen to have a stationary point at x = 1. However, the 
Lagrangian cost is a maximum at the constrained minimum. 
Next, the penalty function is added to the Lagrangian cost. 
In this case, we let P = 5. Now the augmented Lagrangian cost 
has a stationary point at the constrained minimum and the 
second derivative is positive. Note that all three costs have 
the same value at the constrained minimum. 


43 SIGNIFICANCE OF THE LAGRANGE MULTIPLIERS 

Not only are Lagrange multipliers an effective and 
elegant method of including constraints into the optimization 
problem, but as with many problems solved with Lagrange 
multipliers, they provide special information about the system 
being optimized. To show this, consider the optimization 
problem with equality constraints only. At the constrained 


minimum 


vj - + )T 7 2 = ° 


(4.21) 


Hence, it must be that 


^ „yr 7f (4.22) 

But the change in cost for a small change in x is 

A <j> = A 2 (4.23) 


Substitution of Equation (4.22) into (4.23) gives the result 
that 



(4.24) 


Therefore, it is seen from Equation (4.24) that the Lagrange 
multiplier indicates the sensitivity of cost to a change in 
constraints. A similar result holds for the inequality 
constraints . 


4.4 UNCONSTRAINED MINIMIZATIONS 

Optimization problems can be divided into two distinct 
categories : constrained function minimizations and 
unconstrained function minimizations. In Section 4.6, the 


method of mathematical programming via augmented Lagrangians 
£10] will be discussed. With this technique, the constrained 
optimization problem is solved as a series of unconstrained 
optimization problems. In this section, and in Section 4.5, 
the methods used for solving the unconstrained portion of the 
constrained optimization problem are discussed. 

Suppose one wishes to minimize the cost function <J> not 
subject to any constraints. The goal is to find that x which 
produces a minimum <J> . A necessary condition for function 
minimization is that the gradient of the function be equal to 
zero- The simplest method of searching for a minimum <P in 
x-space is to first start at some point and evaluate the 
gradient at this point. Since we are looking for a minimum, 
and hence a decrease in (J>, the negative of the gradient is 
taken as the search direction. Then along this semi-infinite 
line, a minimum will exist. A line search is carried out to 
find the location x^ + ^ of this minimum. Then the gradient of 
g K +iis found at this new point and the entire process is 
repeated until V<f> is approximately equal to zero. This is 
known as a steepest descent gradient search. 

Although the steepest descent search is very simple to 
implement, the convergence to the minimum can be very alow, 
especially when used in conjunction with the method of 
mathematical programming via augmented Lgarangians. Hence, a 
more elegant search routine is needed . The method chosen for 
this research is attributed to Broyden [20] and is one of a 
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class of variable metric methods/ also know as quasi-Newton 
methods. The motivation for these methods is illustrated by 
the Taylor series expansion of the gradient of the cost. 

4 3 K (4.25) 


where 4 Z(< . 

The vector operators in Equation (4.25) are 



(4.26) 


(4.27) 


The matrix VV0 is the so-called Hessian matrix, A. To 
determine the location where the gradient of the cost is zero 
in this Taylor series approximation implies taking a step Ax 
where 



(4.28) 


Notice that this is just a Newton-Raphson step. Figure 4.3 
shows graphically the benefits of using quasi-Newton procedure 
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near the minimum. First of all, the gradient is the vector 
which points in the direction of greatest slope. It does not 
point in the direction of the minimum. Secondly, the gradient 
contains no information about the step size one should take to 
get to the minimum along the gradient search direction. On 
the other hand, near a minimum, a Newton-Raphson search 
direction vector points directly to the minimum, both in 
direction and magnitude. Hence, a search which uses & 
variable metric procedure is usually much more efficient, 
albeit more complicated, than a simple gradient search. 

There exist many schemes for determining the inverse 
Hessian matrix H by iteration. These schemes use the position 

'V 

vector and the gradients at current and previous steps to 
iterate on the H matrix. The method used in this study is 

'Xi 

attributed to Broyden C20] and is one in a general class of 
variable metric methods formulated by Broyden. To perform a 
minimization of 4> using Broyden 1 s method, one first chooses 
some positive definite matrix as an initial guess. Often, 

'Xi 

for convenience, this is a multiple of the identity matrix. 
Next the search direction d is found. 

cL * - H k (4.29) 

This direction and the point x v define the line 

A A 

% = X K + 


c * o 


(4.30) 


where c is the scalar step size. Along this semi-infinite 
line, a minimum is found, or at least a significant 
improvement in <j> . This new point is . The gradient at 


this position is g, 


K+l 


The new H matrix is then taken to be 


H 


K+l ~ W* * 


( 


i + a * t « 

a* 


) 


__ A A<^ Hk 


(4.31) 
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AX A 


where 


This equation is known as the Broyden-Pletcher-Shanno (BFS) 
formula. If the cost function <f> is quadratic, it can be 
shown that H and the minimum will be found exactly in n or 
fewer steps if exact line minimizations are performed during 
each line search. It is assumed that these good qualities 
will apply to a non-quadratic function sufficiently near a 
minimum . 

There are several useful qualities of the BFS variable 
metric search routine. First, only the value of <P and its 
first derivatives need to be found. This saves both 
analytical work and computation time since second derivatives 
are not explicitly evaluated as in a true Newton -Raphs on 
search. Secondly, unlike other variable metric methods, 
notably the method of Davidon, Fletcher, and Powell (DFP) 
[11], exact line minimizations need not be found. In the DFP 
method, the H matrix may become semi-definite near a minimum 


if line minimizations are not performed very precisely. To 
find the exact minimum along a line can be computationally 
expensive and, hence, the BFS routine was chosen since line 
minimizations need not be computed exactly to insure a 

positive definite H matrix and a finite convergence rate [21]. 

a. 

4.5 STEPS IZE CONTROL 

At every step of the search routine a search direction 
is determined. Then, along this line, the vector x is a 
function of the step size c as given by Equation (4.30). Along 
this line, <P may be considered a function of the single 
scalar c. F-.t/.are 4.3 shows the contours of a sample cost 
function in x-spaee. The vector starting at x^ is the search 
vector. One must search, that is vary c through a dicrete set 
of values, to find a step size which produces an acceptable 
reduction in <f>(c). Naturally, one would like to determine 
that step size with the least amount of computational effort. 
A line search for a local minimum of <J> can be very time 
consuming since <{> must be evaluated many times to find the 
minimum precisely. 

Dixon [213 has shown that exact line minimizations may be 
unnecessary to achieve good convergence times when using the 
BFS method. In fact, some of the best results (i.e., fastest 
convergence to the minimum) were obtained while using an 
acceptable step size rule. The Armijo rule [22] used in this 
investigation is nearly identical to the rule used in Dixon’s 
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study . 

The Armijo rule says that step size used at each 
iteration will be c such that 

— V> m 

C* c.i (4.32) 


where y ia a positive number less than unity and c satisfies 
the conditions that 

~ < ~ (4.33) 

Z > o 

and ra is the smallest nonnegative integer which satisfies the 
condition that 


<!>(%*) - -<rc y , *7^(20 T ifc 


(4.34) 


For this study, Y * 0*35 and cr a 0.45. 

Figure 4.4 shows graphically the interpretation of the 
Armijo rule. The curved line is the left hand side of 
Equation (4.34) while the straight line is the right hand 
side. For c to be an acceptable step size, the difference 
between <f>(c) and <<>(0) at the point c must be less than the 
value of a line drawn through the origin with a slope equal to 
a times 3<f>(0)/3c . This can always be satisfied if the step 
size is made small enough. One can see clearly that although 
this rule does not attempt to find an exact line minimum, it 


does force the search routine to move closer and closer to a 
minimum with each succesive iteration since only improvements 
in the cost function are accepted. 

In summary, the method of solving for the unconstrained 
minimum has tv\> main parts : the search direction phase and the 
line search phase. At each step of the search, a search 
direction is found by preraultiplying the negative of the 
gradient with the current estimate of the inverse Hessian. 
Once the search direction has been determined, the Armijo rule 
is used to find the stepsize in that direction which provides 
a significant reduction in the cost function. This procedure 
is repeated until the gradienv of the cost is zero, indicating 
a minimum has been found. 


4 6 CONSTRAINED OPTIMIZATIONS 

As mentioned earlier, the constrained optimization 
problem can be solved as a series of unconstrained problems. 
The general procedure is to form the augmented Lagrangian cost 
function by adding penalty functions to the Lagrangian cost. 
These penalities, in essence, add a large cost to the 
Lagrangian cost if the constraints are not met. This tends to 
force the solution to satisfy the constraints. 

The augmented Lagrangian cost is 






* fl *■ fs 
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(4.35) 
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where i a indicates active inequality constraints 

i^ indicates inactive inequality constraints 
P is a positive number 
W is a positive number 

Hence the gradient cf the augmented Lagrangian cost is 

7?!' Ay * l T H *PT r vr (4 36) 

+ V£ S; 7 S; + w£(»;-|9,|)78; 

i*ik 

The procedure for solving for the constrained minimum ia 
as follows. Initial guesses are chosen for the Lagrange 

multipliers v and n . Initial values for the penalty 

a, 

function coefficients P and W are also selected. The 
augmented Lagrangian cost function is then considered an 
unconstrained function of x only . This function is then 
minimized by an appropriate unconstrained minimization 
procedure. What is meant by appropriate will be discussed 
shortly. Once the minimum of J is found for the given values 

a 

of v , n , P, and W, the Lagrange multipliers and penalty 
function coefficients are then updated and the entire process 
is repeated until convergence. 

The iterative updating of the Lagrange multipliers occurs 
immediately . after each unconstrained minimum of the augmented 
Lagrangian cost. To demonstrata how this procedure works, 
consider the gradient of the Lagrangian and augmented 
Lagrangian costs with equality constraints only. The gradient 
of the Lagrangian cost is 
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The gradient of the augmented Lagrangian cost is 


vji - 7 ^ + )/ T vt + pr r vr 
, * (?+ ?/>/ 


(4.38) 


In the limit of infinite P, with the Lagrange multipliers set 
to zero, the minimum of the augmented Lagrangian cost is a 
constrained minimum of the objective cost function. The value 
of Piji acts like a Lagrange multiplier and in fact will equal 
the true multiplier in the limit of infinite P. For large but 
finite P, we expect that the quantity v + P^ will be 

'V. 

approximately equal to the true Lagrange multiplier at the 
minimum of J . Hence, the update for the equality Lagrange 
multipliers is 


. p? 


{4.39) 


Similarly, the update for the inequality constraint 
multipliers is 

( V(9i - I I) if <o 
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(4.40) 
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Immediately after each update stage/ the Kuhn-Tucker 
conditions for optimality are checked. If these conditions 
are satisfied, then the iteration procedure is stopped since 
the necessary conditions for optimality have been met. Of 
course, this does not guarantee that the constrained 
stationary point will be a minimum. But in practice, this is 
generally the case. 

Note that to solve the constrained minimization problem, 
one need not work with Lagrange multipliers at all. One can 
simply pick P and W to be very large. Then the minimum of the 
cost plus the penalty function will be approximately equal to 
the constrained minimum of the cost function. This method of 
solution, however, will not work well due to the extremely 
slow convergence rates one encounters in the search for a 
minimum. In a two-dimensional problem, one can imagine that 
the augmented cost function is an elevation map. The valleys 
created by the use of penalty functions can be very narrow 
with steep walls. Gradient methods tend to search across the 
valleys and not along them. Hence, convergence to a minimum 
can be very slow. Even the variable metric methods will not 
work well except for regions very close to the minimum. 

However, if the Lagrange multipliers are known 
approximately, the penalty functions need not be so severe and 
the variable metric methods tend to converge faster. This is 
the motivation for the technique of nonlinear programming via 
augmented Lagrangians. Note that a steepest descent search is 
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still not appropriate, however, because the augmented cost 
still will show some of the poor conditioning of contours as 
dicussed above. A variable metric method is required for good 
convergence . 

The values of the penalty functions can be increased with 
each iteration to further increase the requirement that the 
constraints be met. By starting the values of W and P at 
moderate levels at the start of the optimization and 
increasing them at each iteration, one takes full advantage of 
the method of augmented Lagrangians . A simple scheme for the 
update for P and W is given in [11]: 



V 


kP M 


(4.41) 


(4.42) 


where k is a number greater than or equal to one. 

In summary, the rotor optimization problem was formulated 
and the method of solution was outlined. A cost function was 
defined which is a measure of the severity of mistune in the 
rotor. The stability requirements of the problem were 
introduced as constraints which must be satisfied. The cost 
function is to be minimized subject to these contraints. 
Finally, a review of the optimization theory used in this 
study was presented. The results of this optimization 
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procedure will be given in Chapter 5, 


5 . RESULTS OF THE CONSTRAINED OPTIMIZATION PROCEDURE 

In this chapter, the optimization procedure described in 
Chapter 4 is used to determine optimal mass mistune patterns. 
When implemented, these patterns provide the greatest 
stability margin at the aeroelastic design point for the least 


amount of 

mistuning. 

In 

Section 5.1, 

the 

procedure 

for 

optimizing 

the rotor 

is 

outlined and 

the 

behavior of 

the 


optimized rotor at its aeroelastic design point is presented. 
In section 5.2,- the off-design behavior of the rotor is 
analyzed. The next two sections address the issue of the 
actual implementation of mistuning in a rotor. Although the 
designer may specify a certain mistune pattern, the actual 
mistune pattern of the rotor will be different due to 
manufacturing tolerances and changes in the natural 
frequencies of the blade which occur over the life of the 
blade, due to its operating environment. In Section. 5.3, the 
sensitivities of the stability margin to small errors in 
mistune are examined. Finally, the use of optimal mistune 
patterns as a guide to design of near-optimal rotors is 
considered in Section 5.4. For a 14-bladed rotor, the 
optimally mistuned rotor will in general have 14 different 
natural frequencies of the 14 blades- In this section, the 
optimal mistune pattern is approximated by two, three, and 
four tones or frequencies of blades. 


5.1 OPTIMAL MISTUNE PATTERNS FOR A ROTOR AT ITS AEROELASTIC 
DESIGN POINT 

The rotor which was optimized in this study is the same 
rotor which was first introduced in Section 2.5. The 
geometric and aerodynamic properties at the typical (35 
percent span) section are given in Table 2.1. 

The procedure for optimizing the rotor in this 
investigation was as follows. First, an initial guess was 
chosen for the optimal mass raistune pattern. -Also, initial 
guesses were chosen for the Lagrange multipliers. A 
constrained minimum was then found for the case requiring a 
minimum damping ratio greater than or equal to C = -0.C05. A 
second set of constraints requires that all mass changes be 
positive. It should be noted that there are in fact many 
constrained local minima. Which minimum found depends on many 
factors including the initial guess of the optimal mistune, 
the initial choice of Lagrange multipliers, t) , the values of 
the penalty function coefficients, W, the details of the step 
size rule, and the method used to perform the unconstrained 
optimization phase of the algorithm. Only one local minimum 
was found in the search for optimal mistune patterns of the 
14-bladed rotor studied in this report. No attempt was made 
to find other local minima. However, a 13-bladed case and a 
12-bladed case were also examined. In the ] 2-bladed case, 
several local minima were discovered. That work will not be 
presented in this report except to say that the minima found 
in these cases had approximately same level of misnune as the 


optimally mistimed 14-bladed case. It does not appear, 
therefore, that a rotor with an odd number of blades performs 
significantly better or worse than a rotor with an even number 
of blades when the rotors are optimally mistuned. 

As the next step, the optimal mistune pattern for a 
minimum damping ratio greater than £ = -0.004 was found. The 
initial guesses for the optimal mass mistune and the Lagrange 
multipliers were taken to be the result of the optimization of 
the £ = -0.005 optimization problem- The optimal mistune 
patterns for the cases of £ = -0*003, -0.002, -0.002, 0.0, 
0.001, 0.002 were found sequentially in a completely analogous 
fashion, using the preceding optimization results for the 
initial conditions of the next successive optimization. Of 
these cases, the \ =* 0.002 case did not fully converge. The 
Lagrange multipliers at this point were not found, and the 
actual damping ratio produced by the partially converged 
solution was £ =* 0.00188. 

It was not possible to obtain converged solutions to the 
constrained optimization problem for stability margin 
requirements greater than £ = 0-001- The difficulty arises in 
the evaluation of the eigenvalues by integration of the 
derivatives of the eigenvalues and eigenvectors. As the 
system is mistuned, some of the eigenvalues of the system may 
become very close to one another. This is especially true of 
those eigenvalues which lie on the stability margin 
constraint. When the eigenvalues become very closely spaced. 


the first order approximation to the eigenvalue derivatives is 
valid in a very small region of mass mistuning space and the 
accuracy of the integration becomes poor. Furthermore, if two 
eigenvalues should become equal or very nearly equal, the 
integration will break down altogether. The integration 
scheme was not sophisticated enough to avoid these pitfalls by 
dynamically adjusting the integration step size or to avoid 
singularities by rerouting the integration path around them. 
This turned out to be the limiting factor on the usefulness of 
the optimization procedure. 

As stated in Chapter 4, the goal of this analysis is to 
determine the optimal mistune pattern which produces a 
required minimum damping ratio. Figure 5.1 shows the mistune 
patterns found for the eight cases described above: 
5- -0.005, -0.004, ..., 0.001, 0.00188. 

Consider the case of 5 =* -0.005 (see Figure 5.1). This 
type of mistune pattern is called almost alternate mistuning. 
Notice that the odd numbered blades have no change in mass 
from their nominal mass (i.e., E ^ ss 0 ) • Tlie ev en blades all 
have nearly nearly equal masses except for blade 2, which has 
a mass mistune of zero, and blade 14, which has a mass mistune 
of about half that of the other even numbered blades. 

This almost alternate mistune pattern grows in magnitude 
for increasing ? but does not change in nature until 
C = 0.0. At this value of required damping ratio, the mass of 
blade 9 becomes nonzero. At C =» 0.001, the mass of blade 



number 2 becomes nonzero but still very small compared to the 
mass of the other even numbered blades . 

Notice that the mistune patterns shown in Figure 5.1 all 
resemble, to some degree, the truly alternate mistune pattern 
where all the odd numbered blades have 0 and all the even 
numbered blades have equal nonzero e^. One might expect, 
therefore, that although the mistune patterns found in this 
investigation may be optimal, they are not ^ nificantly more 
cost effective than the truly alternate mistune pattern. 
However, as demonstrated by Figure 5.2, this is clearly not 
the case. The figure shows the cost of mistuning (as defined 
by Equation (4.3)) versus the stability margin achieved by the 
mistuning. The upper curve is the cost of the truly alternate 
mistune pattern and the lower curve is the cost of the optimal 
mistune pattern. The optimal mistune pattern is seen to 
deliver much more stability for a given level of mistune or, 
alternatively, a much lower cost for the same level of 
stability. To achieve a damping ratio of 0.00188, the optimal 
mistune pattern requires about 45 percent less mistuning than 
alternate mistuning. 

Some insight into why the optimal mistune patterns are so 
effective can be gained by looking at the eigenvalues in the 
complex plane. Figures 5.3a-i show the eigenvalue for 
C = -0.00602 (tuned), -0.005, -0.004, -0.003, 0.001, 
0.00188. As the stability margin becomes greater, more and 
more eigenvalues just barely satisfy the constraint that all 
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eigenvalues lie to the left of the ray emanating from the 
origin with a damping ratio of C . For a damping ratio of 
C « 0.00188 (see Figure 5.3i), four of the 14 eigenvalues lie 
on the stability margin constraint. This indicates that the 
mistune pattern is very efficient, since it does not do 
unnecessary work by pushing some of these eigenvalues further 
to the left than required. In contrast, only one of the 14 
eigenvalues of the alternate mistune pattern will, in general, 
lie on the constraint. 

As a by-product of the optimization procedure, the 
Lagrange multipliers of the active constraints are determined. 
Recall from Section 4.3 that the Lagrange multipliers indicate 
the change in the optimal cost for a small change in the 
constraints. Hence, by summing the negative of the Lagrange 
multipliers associated with the damping ratio constraints, the 
local slope of the cost versus damping ratio curve (see Figure 
5.2) is determined. To check this result, the slope of the 
optimal cost curve was determined by using second order finite 
difference operators. These results are plotted in Fi.gure 
5.4, along with the slope predicted by the Lagrange 
multipliers. Note the generally good agreement between the 
two. The difference between the two can be attributed to the 
relatively large A? used in the differencing and the limited 
accuracy to which the Lagrange multipliers were computed. 


5.2 OFF DESIGN PEFORMANCE OF OPTIMALLY MISTUNED ROTOR 

The rotor has been mistimed to achieve given stability 
requirements at the aeroela&tic design point of the rotor. It 
remains to be seen if this mistuned rotor will be stable over 
the entire operating range. To determine this, a modified V-g 
diagram was constructed. The negative of the damping ratio of 
the rotor is plotted versus the reduced velocity. For 
stability, the damping ratio must be positive, i.e., the curve 
must lie below 5 = 0. Two cases were examined. In the first 
case, the reduced velocity was varied while all other 
parameters, including the Mach number, were held constant. In 
the second case, the Mach number and reduced velocity were 
varied together to simulate a fan rotor running up its 
operating line. 

Plotted in Figure 5.5 are three important stability 
curves for the constant Mach number case. The upper curve is 
the tuned damping ratio plotted against the relative reduced 
velocity of the rotor. Bendiksen E 3 3 has shown that for 
rotors which can be modelled well with only one degree of 
freedom per blade, the tuned position is always the least 
stable. The lower curve is the damping ratio of the the 
centroid of the tuned eigenvalues. Recall from Section 2.4 
that this is to first order the best damping ratio that one 
can achieve by mistuning. Hence, the mistuned damping ratios 
should lie between th^Le two curves. Finally, the center 
curve is the rotor which has been optimally mistuned at the 
aeroelastic design point. As expected, this curve lies 


between the tuned deunping ratio curve and the blade self 
damping curve. 

As Figure 5 . 5> clearly shows, the rotor has been 
stabilized for all reduced velocities less than or equal to 
the aeroelastic design reduced velocity. For the worst case 
(the tuned case), the rotor flutters at a reduced velocity of 
V * 1.7. The best case (the centroid of the eigenvalues) 
produces a flutter speed of V = 2.9. The use of optimal 
mistuning has increased the flutter speed from the worst case 
to V =» 2-05. Hence, the flutter speed has been increased by 
about 20 percent over the tuned flutter speed. In the best 
possible case, one could mistune the rotor to achieve a 70 
percent increase in flutter speed. 

In this example, the flutter speed of the rotor has been 
increased by 20 percent. This increase in flutter speed can 
prevent a rotor which would flutter in its tuned state from 
fluttering at its aerodynamic design point. On the other 
hand, if a rotor does not flutter in its tuned state, the 
increase in flutter speed could be used to reduce the chord of 
the blades without inducing flutter due to the subsequent 
increase in reduced velocity. 

As a second example the case of a given fan running up 
an operating line was considered. In this case, the Mach 
number and the reduced velocity were held proportional to one 
another. Figure 5.6 shows a "V-M-g diagram" (damping ratio 
versus reduced velocity and Mach number) . Again the upper and 
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lower curves are the tuned damping ratios, and the damping 
ratio of the centroid of the tuned eigenvalues, respectively. 
The middle curve is the optimally mistuned rotor. Notice that 
again, the mistuned damping ratios lie between the tuned 
damping ratio and the centroid damping ratio. In this 
example, however, although the rotor was stabilized at a 
reduced velocity of 2.02 and a Mach number of 1.317, the rotor 
appears unstable at reduced velocities and Mach numbers lower 
than the aeroelastic design point. Note that there is a 
single region of instability just before the operating point 
is reached. It is possible that if the rotor had been more 
severly mistuned, the mistuned damping ratio curve would not 
have gone from positive to negative before the operating point 
is reached. Unfortunately, as previously dicussed, the 
optimization procedure failed to converge for damping ratios 
greater than C * 0.001. 

5.3 SENSITIVITY TO ERRORS IN MISTUNING 

An important issue, which must be addressed before an 
optimally mistuned rotor is actually used in aeroelastic 
experiments, is the question of sensitivity to manufacturing 
errors. Although the designer may specify a certain mis tune 
pattern, he must accept the fact that in the manufacturing 
process there will be certain tolerances which cannot be 
obtained. Hence, the actual mistune pattern which is 
implemented will be somewhat less than optimal. 


The actual 


mistime pattern will be 
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where e. is the error in mistuning the rotor. 
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To investigate this problem, errors in mistuning were 
introduced into the optimally mistuned rotor with a stability 
margin of 0.00183. The gradient of the damping ratio of each 
eigenvalue was found with respect to the mistune variables 

e . . The actual mistune patterns were then taken to be 
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In Equation (5.2), the term — vc n /|J VC n j| L 8 the vector which 
points in the most destabilizing direction of the nth 

eigenvalue. This vector is normalized to have a length of 
unity. The vector is multiplied by E/fT, where E is the root 
mean square of the entries of the mistune error vector, then 
added to the nominal mistune pattern. A typical value of E 
for an actual rotor is about 0.01. The rotor was mistuned 
using the pattern specified by Equation (5.2) with the N 
different gradients corresponding to the N different 
eigenvalues. Then the case with the worst damping ratio after 
the errors had been introduced is, to first order,- the worst 
possible case for a given value of E.. It was found that for 
E =* 0.01, the stability is reduced from 0.00188 to -0.00317 
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(see Figures 5.7). Hence it is Been that the optimal mistune 
pattern for 5 *» 0.00188 is very sensitive to small changes in 
raistuning . 

The eigenvalues of the perfectly mistuned and the worst 
case imperfectly mistuned systems are plotted in Figure 5.8. 
In this case, the eigenvalue which moves the most to the right 
in the complex plane is the eigenvalue on the root locus of 
the tuned 77.14° interblade phase angle eigenvalue . Before 
the introduction of mistuning error, this is one of the 
eigenvalues which lies on the stability margin constraint. 

In a sense, the optimal mistune pattern is very sensitive 
to errors because it is an optimal mistune pattern. Figure 
5.7 shows that the optimal cost curve has a very shallow slope 
at f =» 0.00188. This implies that for a small increase in 
raistuning, a large improvement can be made in the stability 
margin. But for this same reason, a small change in mistuning 
can graatly reduce the stability margin. 

This same sort of sensitivity analysis was carried out on 
the alternately mistuned rotor with a perfectly mistuned 
stability margin of ? ® 0.00171. Again the error vector was 
chosen to be in the worst possible direction for an error 
vector with a root mean square of 0.01. The stability margin 
in this case was degraded from 0.00171 to 0.00047 as shown in 
Figure 5.7. As it turns out, alternate mistuning is less 
sensitive to mistuning errors than optimal mistune. By the 


symmetry of the problem, the gradients of the damping ratios 
all are in the direction of an alternate mistiming direction. 
Since alternate mistuning does not give large improvements in 
the minimum damping ratio for a small change in mis tune, one 
would not expect the errors in mistuning to produce a large 
change in the stability margin since the worst errors are in 
an alternate mistune direction. 

An interesing result of the perturbation analysis 
presented in Chapter 3 is that small amounts of mistuning do 
not significantly change the stability of the system from the 
originally tuned configuration. To show this, recall that the 
derivatives of the eigenvalues of the system are: 
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For the case of mass mistuning. Equation (5-3) becomes 
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It can be shown that, for the tuned case, the quantity on the 
right hand side of Equation (5.4) is equal to -1/N. Hence, 
for small mass perturbations about the tuned position, the 
change in an eigenvalue is equal to 
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Therefore, k >ie eigenvalues of the slightly mistuned system 


depend only on the net mass added to the system, and not on 
the mistune pattern- But if equal masses are added to all the 
blades , one would not expect the damping ratios of the 
eigenvalues to change significantly since the system is still 
tuned. Therefore, the stability margin of the tuned system is 
insensitive to small amounts of mistune, no matter what 
mistune pattern is used. This is shown clearly in Figure 5.7. 
Both the optimal mistuning and the alternate mistuning cost 
curves are very steep at the tuned postion. Relatively large 
amounts of mistune are required to cause small changes in 
stability about the initially tuned position. 

5-4 TWO, THREE, AND FOUR TONE APPROXIMATIONS TO OPTIMAL 
MISTUNE PATTERNS 

The cost function minimized in this optimization reflects 
the practical difficulty associated with reaching a certain 
magnitude of mistuning. However, the difficulty with 
constructing a mistuned rotor is not only in the level of 
mistuning, but also in the complexity of the mistune pattern. 
For example, for the rotor examined in this study, the optimal 
mistune pattern required that there be nine different blade 
natural frequencies for the 14 blades on the rotor. Hence, to 
actually construct an optimally mistuned rotor, one would have 
to construct approximately N/2 different types of blade for a 
single rotor. This would prove to be very costly. In this 
section, several suboptimal mistuning patterns are created by 
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approximating the optimal miatune patterns by two, three, and 
four tone miatune patterns, and the performance of these 
miatune patterns are presented. 

The procedure for picking the approximations to the 
optimal miatune was to use the \ =* 0.00188 optimal pattern as 
a guide to provide insight for picking suboptimal miatune 
patterns. First, the two tone approximations were chosen. In 
Table 5.1, the optimal miatune pattern and the two tone 
approximations which were - investigated are given. Of course, 
the goal of these suboptimal patterns is the same as the goal 
for the optimal patterns, i.e., to provide the greatest 
stability margin for the lowest cost. Figure 5.9 shows the 
cost versus stability margin of the two tone patterns. Note 
that two of the patterns result in costs which lie between the 
optimal and alternate mistune costs. The third is about the 
same cost as the alternate mistune pattern. 

Next, several three tone approximations were examined. 
These patterns are presented in Table 5.2. All of the mistune 
patterns tried had a lower cost per stability margin than the 
alternate mistune pattern as seen in Figure 5.10, although not 
dramatically lower. This is an indication that the fine 
detail of the optimal mistune pattern is important. Large 
amounts of mistune will not be effective in preventing flutter 
if this detail is missing. 

Finally, four tone approximations were examined. Table 
5.3 lists the four tone patterns used to approximate the 


optimal mistune pattern. In Figure 5.11, it is seen that 
again, the four tone patterns perform slightly better than 
alternate mistuning, but not as well as optimal raistuning. 
Even with four tone approximations, there is not sufficient 
similarity to the optimal mistune pattern to achieve results 
that are nearly optimal. 

In this chapter, it has been shown that optimal mistune 
patterns can achieve a given stability margin for a relatively 
low level of mistuning . Approximations to these optimal 
mistune patterns using two, three, and four discrete blade 
frequencies, however, perform only slightly better than 
alternate raistuning. Furthermore, optimal mistuning is not 
very robust to small mistuning errors. Alternate mistuning, 
on the other hand is relatively insensitive to errors in 
mistuning. For these reasons, it appears that a practical 
mistune pattern for implementation is the alternate mistune 
pattern . 


6 . CONCLUSIONS 


1. The aeroelastic poles or eigenfrequenci.es can be thought 
of as forming a pattern around their centroid. The location 
of the centroid is controlled by the average blade mechanical 
properties and the aerodynamic blade self damping. It was 
shown that a necessary but not sufficient condition for 
aeroelastic stability is that the blades be self damped. 

2. The distribution of the poles about the centroid is due to 
the unsteady aeroelastic influence of the neighboring blades. 
This pattern of poles can be modified by mistuning the rotor, 
increasing the stability of the less stable poles. However, 
mistuning does not introduce additional damping into the 
system, since the damping ratio of the centroid is unaffected. 
Mistuning makes use of the existing damping to stabilize the 
Litor by decreasing the blade to blade aerodynamic influences, 
cmiteby increasing the stability of the least stable poles. 

3. There are two main mechanisms which can lead to the onset 
of flutter in transonic fans: The first is the lost of the 
blade self damping, as in the case of high incidence stall 
flutter. In such instances, the centroid and the entire 
pattern of poles shifts to the right in the complex plane 
eventually causing the least stable poles to become unstable. 
The second mechanism is due to the increasing destabilizing 
effect of the neighboring blades with increasing reduced 
velocity. The off diagonal influence coefficients, which 


reflect the influence of one blade in the cascade on its 
neighbors , cause the eigenvalues to spread out away from the 
centroid of the poles. Some of the eigenvalues will become 
more stable than the centroid, and some will become less 
stable than the centroid. Hence, in the presence of cascade 
effects, the least stable eigenvalue will necessarily be less 
stable than the centroid of the eigenvalues, the value of 
which is determined only by the blade self damping terms of 
the influence coefficient matrix. 

4. The unsteady aerodynamic forces, derived in terms of 
travelling wave coordinates, can be transformed into a form 
which expresses the force on each blade explicitly in terms of 
the motion of the other blades in the cascade. This linear 
transformation is simply a Fourier decomposition of the forces 
as expressed in the travelling wave coordinates. The use of 
this transformation on both analytical cascade models and 
experimentally measure unsteady aerodynamic coefficients 
generally reveals that the dominant forces acting on a blade 
arise from the motion of the blade itself and its two adjacent 
neighbors. Hence, one would expect that any effective 
mistuning scheme will minimize the influence one blade has on 
its neighbors. This heuristic argument suggests that the 
alternate mistune pattern should be effective. 

5. In an effort to better understand the mechanisms of 
m.istuning, an inverse design procedure was developed which 
determines the optimal pattern of blade structural mistuning 
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for a required increase in aeroelastic stability margin. The 
optimal mistune pattern can achieve a given stability margin 
with a significantly lower level of mistuning than with the 
simple alternate mistuning. The success of optimal mistuning 
does not depend on an even number of blades being present, as 
similar results are found for rotors with odd and prime 
numbers of blades. 

6. The optimal mistune pattern appears to have three salient 
features: First, in all optimal mistune patterns there is 
seen to be some features of the alternate mistune pattern, 
i.e., nearly every other blade is mistuned. This componant of 
the optimal mistune pattern serves to disrupt the dominant 
aerodynamic effects of the neighboring blades. Second, there 
are "break points" around the rotor which disrupt the 
alternate mistune pattern. It is thought that these break 
points prevent longer wavelength disturbances from travelling 
around the rotor. Third, those blades that are mistuned do 
not all have exactly the same amount of mistune. Rather, 
there is a subtle structure to the mistuning which is not 
possible to predict a priori. 

7 . The subtle detail in the mistune pattern appears to be 
very important to the effectiveness of the optimal mistuning. 
The optimal mistuning patterns were found to be very sensitive 
to small errors in mistuning due to the loss of this detailed 
structure. Alternate mistuning, on the other hand, was found 
to be relatively insensitive to errors in mistuning. 
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8. Another consequence of the importance of the subtle detail 
in the optimal mistune pattern is that any practical 
implementation of the mistune pattern using only a small 
number of different blade frequencies will not faithfully 
reproduce all the important features of the optimal mistune 
pattern. It was found that the practical suboptimal mistune 
patterns do not perform significantly better than alternate 
mistuning. 

9. Three distinct regions of mistuning influence were 
identified. Starting from the tuned configuration, there is a 
first region were the stability boundary is insensitive to the 
addition of mistuning. It is thought that most present stages 
operate in this initial insensitive region, which explains why 
they behave similarly and can be analyzed as tuned rotors. 
After several percent of mistuning has been introduced, a 
region of approximately linear increase in stability with 
increasing optimal mistune is entered. Finally, an asymptotic 
limit on the ability of mistuning to increase stability is 
reached. It is apparent that to provide uniformity of 
performance of rotors in service, rotors should be designed 
assuming small amounts of mistuning (i.e. in the initial 
insensitive region) , or if the rotor is to be deliberately 
mistuned for stability, a large amount of mistuning should be 
introduced so that the rotor operates in the latter region of 
insensitivity . 


10. When studying the off design performance of a mistuned 
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rotor, it is useful to construct diagrams similar to the 
traditional V-g diagrams. The stability of the mistuned 
rotor, at least if the rotor can be modelled with a single 
degree of freedom per blade, will lie between two limiting 
curves. The worst possible case is the case of the tuned 
rotor. The mistuned rotor cannot be any less stable than 
this. On the other hand, the most stability one can achieve 
through the use of mistuning is limited by the blade self 
damping. This is the fundamental limitation to the usefulness 
of mistuning. 

11. In order to perform the optimization discussed in Chapter 
4, it was first necessary to develop a method of evaluating 
the eigenvalue and eigenvectors of the equations of motion 
which retains the identity of each eigenmode in a root locus 
sense. This was done by first determining the derivatives of 
the eigenvalues and eigenvectors as a function of the mistune 
of each blade. These derivatives were then integrated to 
determine the eigenvalues and eigenvectors of a given mistune 
pattern. This method of evaluating the eigenvalues and 
eigenvectors is reasonably efficient when used in the 
optimization procedure. Furthermore, the identities of the 
eigenmodes are not lost in their evaluation as the system is 
mistuned. The principle drawback of the method is that the 
integration scheme breaks down whenever two or more 
eigenvalues become nearly equal. This turned out to be the 
limiting factor in the ability to optimize the rotor. For 


! 


i 

4 

4 

j 


93 


damping ratios of 0.002 and greater, the optimization routine 
did not converge due to problems in evaluating the eigenvalues 
and their derivatives. 
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Table 2.1 Parameters of rotor at typical section at the 
aeroelastic design point. 


Number of Blades 

N 

14 

Solidity 

CT 

1 u 409 

Mach number 

M 

1.317 

Recuced frequency 

k 

0.495 

Location of pitch axis 

a 

0.0 

Mass ratio of reference blade 

u 

181.9 

Radius of gyration of blades 

r 

0.4731 

Stagger angle at typical section 

5 

58.99° 


CT . 
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Table 5.1 Two tone mistune patterns. Patterns 2a, 2b, 

and 2c were picked to approximate optimal mistune pattern. 


Percent mass mistuning of ith blade, 


number 

Optimal 

2a 

2b 

2c 

1 

0.00 % 

0.00 % 

0.00 % 

0.00 % 

2 

0.94 

0.00 

0.00 

6.00 

3 

0.00 

0.00 

0 . 00 

0.00 

4 

4.27 

6.00 

6.00 

6.00 

5 

0.00 

0.00 

0.00 

0.00 

6 

5.81 

6.00 

6.00 

6.00 

7 

0.00 

0.00 

0.00 

0.00 

8 

5.09 

6.00 

'6.00 

6.00 

9 

3.47 

0.00 

6.00 

6.00 

10 

7.05 

6.00 

6.00 

6.00 

11 

0.00 

0.00 

0.00 

0.00 

12 

4.76 

6.00 

6.00 

6.00 

13 

0.00 

0.00 

0 • 00 

0.00 

14 

6.26 

6.00 

6.00 

6.00 
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Table 5.2 Three tone mistune patterns. Patterns 3a, 
3b, 3c, 3d, 3e, 3f, and 3g were picked to approximate optimal 
mistune pattern. 


Percent mass raistuning of ith blade, 


Blade number 

Optimal 

3a 

3b 

3c 

i = 1 

0.00 % 

0.00 % 

0.00 % 

0.00 % 

2 

0.94 

0.00 

0.00 

6.00 

3 

0.00 

0.00 

0.00 

0.00 

4 

4.27 

5.00 

6.00 

6.00 

5 

0.00 

0.00 

0.00 

0.00 

6 

5.81 

6.00 

7.00 

7.00 

7 

0.00 

0.00 

0.00 

0.00 

8 

5.09 

5.00 

6.00 

6.00 

9 

3.47 

5.00 

0 . 00 

0.00 

10 

7.05 

6.00 

7.00 

7.00 

11 

0.00 

0 . 00 

0 . 00 

0.00 

12 

4.76 

5.00 

6.00 

6.00 

13 

0 . 00 

0 . 00 

0.00 

0.00 

14 

6.26 

6.00 

7.00 

7.00 

Blade number 

3d 

3e 

3f 

3g 

i = 1 

0.00 % 

0.00 % 

0.00 % 

0.00 % 

2 

0.00 

0.00 

6.00 

0.00 

3 

0.00 

0.00 

0.00 

0.00 

4 

6.00 

6.00 

6.00 

6.00 

5 

0.00 

0.00 

0.00 

0.00 

6 

7.00 

8.00 

8.00 

9.00 

7 

0.00 

0.00 

0.00 

0.00 

8 

6.00 

6.00 

6.00 

6.00 

9 

6.00 

6.00 

6.00 

6.00 

10 

7.00 

8.00 

8.00 

9.00 

11 

0.00 

0.00 

0.00 

0.00 

12 

6.00 

6.00 

6.00 

6.00 

13 

0.00 

0.00 

0.00 

0 . 00 

14 

7.00 

8.00 

8.00 

9.00 
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Table 5.3 Pour tone mistune patterns. Patterns 4a, 4b, 
4c, and 4d were picked to approximate optimal mistune pattern. 


Percent mass mistuning of ith blade, 


Blade number 

Optimal 

4a 

4b 

4c 

4d 

i » 1 

0.00 % 

0.00 % 

0.00 % 

0.00 % 

0.00 % 

2 

0.94 

2.00 
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Figure 2.2 Geometry of the rotor at the typical section. Note 
the stagger angle and blade numbering definitions. The solidity of 
the rotor, a, is equal to 2b/s. 
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Figure 2.4 Unsteady moment coefficients in influence coefficient 
form. Note the four largest coefficients are L^f and L^, 

Indicating that only neighboring blades have a significant influence on 
a given blade. Note also that L_ - L. , . 
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Figure 2.5 Eigenvalues of the tuned rotor, 
of the 14 eigenvalues are unstable. 
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Figure 2.6 Eigenvalues of the alternately mis tuned rotor. The 
even numbered blades have ■ 0.0 while the odd numbered blades have 
a mass mistune of e^ « 0.1. For this case, the mistuning has stabilized 
an otherwise unstable rotor. 
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Figure 3.1 The path of integration for the evaluation of the 
eigenvalues and eigenvectors is chosen to be the straight line connecting 
the end points, ^ and The path is then subdivided into I intervals. 

Over each interval, a fourth-order Runge-Kutta integration is performed. 
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Figure 3.2 Evaluation of eigenvalues of alternately mlstuned rotor 
via fourth-order Runge-Kutta Integration. Number of Integration steps ** 1. 
Note the very poor agreement between the values determined by integration 
and those determined by EISPACK. 
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Figure 3.3 Evaluation of eigenvalues of alternately mlstuned rotor 
via fourth-order Runge-Kutta Integration. Number of Integration steps - 2, 
Note the very poor agreement between the values determined by Integration 
and those determined by EISPACK. 
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Figure 3.4 Evaluation of eigenvalues of alternately mistuned rotor 
via fourth-order Runge-Rutta integration. Number of Integration steps = 3 
Ten of the 14 eigenvalues determined by integration are in excellent 
agreement with those determined by EISPACK. However, four of the 
eigenvalues are not predicted well with three integration steps. 
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Figure 3.5 Evaluation of eigenvalues of alternately mistuned rotor 
via fourth-order Runge-Kutta integration. Number of Integration steps = 1 
With ten integration steps, the eigenvalues are correctly predicted using 
the Runge-Kutta integration scheme. 
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Figure 4.1 Graphic interpretation of stability margin constraint 
Every eigenvalue must have a damping ratio greater than or equal to 5. 
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Figure 4.2 Relationship between cost, Lagrangian cost, and 
augmented Lagrangian cost. The constrained minimum is x = 1. Note 
that at this point, all three cost functions have the same value. 
Furthermore, the slopes of the Lagrangian cost and the augmented 
Lagrangian cost curves are both zero. The augmented Lagrangian cost 
is positve definite at x « 1. 
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Figure 4.3 Near a minimum, a Newton-Raphson search vector points 
directly to the minimum of the cost function, both in direction and 
magnitude. On the other hand, a steepest descent search vector points 
in the direction of largest change in cost, which is not toward the 
minimum unless the cost contours are circular. Furthermore, the steepest 
descent vector contains no stepsize information. 
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Figure 4.4 Graphical interpretation of the Armijo rule (from 

_ 2 

Reference [22]). The stepsize c “ c Y is chosen since m ■ 2 is the 
smallest positive interger which results in a value of <]i(c) - $(0) which 
is less than . 
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Figure 5.1 Mass mistuning vectors of the optimally distuned rotor. 
Note that both the level of mistuning and the character of mistuning changes 
for increasing stability margin requirements. 
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Figure 5.2 Cost effectiveness of optimal mistuning. Note that 
optimal mistuning can achieve a given stability margin for a low level 
of mistuning compared to the alternate mistune pattern. 
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Figure 5.3b Eigenvalues of optimally mistuned rotor. The required 
stability margin is £ - -0.005. 
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Figure 5.3c Eigenvalues of optimally mistuned rotor. The required 
stability margin is Z “ -0.004. 
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Figure 5.3d Eigenvalues of optimally mistuned rotor. The required 
stability margin is £ ■ -0.003. 
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Figure 5.3e Eigenvalues of optimally mistuned rotor. Ths required 
stability margin is \ * -0.002. 
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Figure 5.3h Eigenvalues of optimally mistuned rotor. The required 
stability margin is C = 0.001. 
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Figure 5.3i Eigenvalues of optimally mistuned rotor. The required 
stability is £ * 0.002. However, for this case, the optimization procedure 
did not completely converge. The stability margin of the partially 
converged solution is £ = 0.00188. 
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Figure 5.4 Slope of optimal cost versus stability margin curve. 

The slope of the optimal cost curve was evaluated from the Lagrange 
multiplies of the stability margin constraints, and by finite differncing 
the optimal cost curve. 
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Figure 5.5 Stability margin of a rotor versua reduced velocity for 
the case o. ■ , onfLant Mach number. Note the optimally mlstuned rotor damping 
ratio lies between the blade self damping and the tuned rotor damping 
ratios. 
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Figure 5.7 Sensitivity of the stability margin of a miatuned rotor 
to errors in mistunlng. The arrows Indicate the loss in stability due to 
mis tuning errors of 1% root mean square introduced in the worst poosible 
direction. Note that optimal mistuning is much more sensitive to these 
errors than alternate mistunlng. 
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Figure 5.8 Sensitivity of eigenvalues of optimally mis tuned rotor 
to errors in mistuning. The errors introduced into the system have a 
root mean square value of 0.01, and are introduced in the direction in 
mlstune space which is most destabilizing. 
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Figure 5.9 Coat effectiveness of two tone mistune patterns. Mistune 
patterns 2a, 2b, and 2c were chosen to approximate optimal mistune pattern 
for l - 0.0018. 


133 


ORIGINAL PAGE 
OF POOR QUALITY 


Cost Effectiveness of 
S uBOPTIMAL MlSTUNING 

O Three tone approximations 
to optimal 

— Optimal miatuning 

Alternate mistuning 


c 

3 ao 
U O 
0 ) « 
•H O 
X 


Stability Margin, £ 


Figure 5.10 Coat effectiveness of three tone mistune patterns. 
Mlstune patterns 3a, 3b, 3c, 3d, 3e, 3f , and 3g were chosen to approximate 
the optimal mistune pattern for Z m 0.0018. 
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Figure 5.11 Cost effectiveness of four tone mistune patterns. 
Mistune patterns 4a, 4b, 4c, and 4d were chosen to approximate optimal 
mistune pattern for C - 0.0018. 


Q 



